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pg I Abstract 

Interference between nodes is a critical impairment in mobile ad hoc networks (MANETs). This 
paper studies the role of multiple antennas in mitigating this interference. Specifically, a network is 
studied in which transmitters opportunistically transmit and zero-forcing beamforming is applied at each 
receiver for canceling the strongest interferers. The performance of this approach is analyzed by evaluating 
the transmission capacity of a network with Poisson distributed transmitters and i.i.d. Rayleigh fading 
channels. This metric corresponds to the maximum density of transmitting nodes subject to an outage 
O , constraint for a given signal-to-interference ratio (SIR). Mathematical tools from stochastic geometry are 

applied to obtain scaling laws for the transmission capacity and characterize the impact of inaccurate 
^ channel state information (CSI). In particular, for small target outage probabilities, the transmission 

(^ ' capacity is proven to increase following a power law, where the exponent is the inverse of the size of 






each node's antenna array or larger, depending on the path-loss exponent. Moreover, CSI inaccuracy 
is shown to only decrease the constant in the scaling law of the transmission capacity as the outage 
[-*.. i probability goes to zero, provided that an appropriate length for the CSI training sequence is used. The 

^^ . needed length of such training is also derived. Numerical results suggest that using merely one additional 

(•^ , antenna at each node increases the transmission capacity by an order of magnitude or more, even when 

the CSI is imperfect. 

X 

o3 ' I. Introduction 

In a mobile ad hoc network (MANET), the mutual interference between nodes poses a fundamental limit 
on the throughput of peer-to-peer communication. One approach for mitigating the effect of interference 
is to provision the nodes in a MANET with multiple antennas and use the spatial degrees of freedom 
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created by these antennas to cancel a portion of the interference. We consider such an approach in this 
paper. Namely, each receiver uses zero-forcing beamforming to attempt to cancel the interference from 
the strongest interferers. Each transmitter simply chooses a random beam and transmits according to a 
threshold rule as in [1]. This approach requires only limited local coordination making it suitable for a 
MANET. A basic question is then to quantify the gains in network performance as a function of the system 
parameters. Furthermore, implementing this approach requires some amount of channel state information 
(CSI) at each receiver. In a MANET environment, there will naturally be inaccuracy in this CSI. A second 
key question is then to determine how this inaccuracy will impact the network performance. 

We provide answers to both of the above questions in terms of the transmission capacity of a simple 
network consisting of Poisson distributed transmitters and spatially i.i.d. Rayleigh fading. This metric, 
introduced in [2], is defined as the maximum density of transmitters so that a typical transmitter will 
satisfy an outage probability constraint for a target signal-to-interference-and-noise ratio (SINR)lll By 
studying the scaling of transmission capacity as the outage probability goes to zero, we are able to 
quantify the gains in network performance due to spatial interference cancellation as a function of the 
relevant system parameters, as well as the effect of CSI inaccuracy. These results suggest that multiple 
antennas can significantly improve the performance of MANETs, even with inaccurate CSI. 

A. Prior Work and Motivation 

As noted previously, our primary performance metric is the transmission capacity. For Poisson dis- 
tributed transmitters, this metric is studied in [2] for single-antenna ad hoc networks assuming fixed 
transmission power and an ALOHA-like medium access control (MAC) layer. Transmission capacity has 
also been used to make a tractable analysis of opportunistic transmissions [1], distributed scheduling [3], 
coverage [4], network irregularity [5], bandwidth partitioning [6], and successive interference cancellation 
(SIC) [7] in ad hoc networks. Additionally, in [1], transmission capacity is analyzed for several models 
of multi-antenna ad hoc networks. The work in [1] differs from this paper in that the antennas are not 
used for interference cancellation, but rather for interference averaging through different multiple antenna 
diversity techniques. 

As discussed in [1], [2], transmission capacity is closely related to the notion of transport capacity 
introduced in [8] and further studied in a number of papers including [9]-[13]. Transport capacity focuses 
on the scaling of a network's total end-to-end throughput per unit distance as a function of the network's 
size, while our focus here is on the number of single -hop transmissions as a function of the transmitter 
density. Furthermore, most work on transport capacity assumes perfect scheduling and zero-outage, while 
here we focus on a random access model with outages. 

'This metric focuses on the one-hop performance of a MANET and does not explicitly account for multi-hop routing. 
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In addition to spatial interference cancellation, there are a number of other possible approaches for 
mitigating interference in MANETs. For example, the interference alignment approach in [14] achieves 
the optimal number of degrees of freedom in a high signal-to-noise ratio (SNR) setting. This approach 
appears daunting in practice because it requires jointly designed precoders and perfect CSI of interference 
channels. In contrast, the approach considered here only requires "local" channel state information and no 
coordination of transmit precoders. Another method for interference management, used in many practical 
MAC protocols, is to create an interferer-free area - a guard zone - around each receiving node through 
carrier sensing. As shown in [3], by optimizing the guard-zone size, this method leads to significant gain 
in a single-antenna MANET's transmission capacity with respect to purely random access. Here, the use 
of interference cancellation, can be viewed as creating an effective guard zone for each receiving node, 
without requiring that other nearby transmitters are suppressed. 

A number of other recent papers have addressed other aspects of multi-antenna MANETs. For example, 
beamforming or directional antennas [15] have been integrated with the MAC protocols for MANETs to 
achieve higher network spatial reuse or energy efficiency [16]-[28]. In addition, multi-antenna techniques 
have been shown to improve the efficiency of routing protocols for MANETs [29]-[33]. Directional 
antennas have been studied for suppressing interference in MANETs by spatial filtering [23]-[27]. 
Directional antennas are only suitable for environments with sparse scattering. In contrast, beamforming 
is suitable for both sparse and rich scattering, and is hence adopted in this paper as well as in [18]- 
[20], [28] for spatial interference cancellation. Most prior work focuses on designing MAC protocols 
and relies on simulations for investigating network throughput [16]-[27]. The capacity of MANETs with 
beamforming or direction antennas are analyzed in [34]-[36]. In [34], [35], the use of directional antennas 
are shown to increase the linear scaling factor of network transport capacity. In [36], the transmission 
capacity for multi-antenna MANETs is analyzed, where interference is treated as noise and suppressed 
by averaging through beamforming. In view of prior work, there still lacks of theoretic characterization 
of the relationship between the transmission capacity of MANETs and spatial interference cancellation. 
Furthermore, the important issue of how CSI inaccuracy affects the throughput of MANETs has not been 
analyzed in [16]-[27], [34]-[36]. 

B. Contributions and Organization 

Our main contributions are summarized as follows. This paper targets a MANET with single-stream data 
links and perfect symbol synchronization between nodes. First, we give an approach for using zero-forcin 
beamforming to cancel interference in a MANET and thereby improve network transmission capacitjt 

^The zero-forcing method is used for analytical simplicity, and the extension to minimum-mean-squared-error (MMSE) 
beamforming is straightforward. 
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The number of canceled interferers depends on the number of antennas at each receiver. Moreover, transmit 
beamforming vectors are randomly selected to avoid iterative receive beamforming and potential network 
instability. Second, assuming Poisson distributed transmitters and spatially i.i.d. Rayleigh fading channels, 
bounds on the signal-to-interference ratio (SIR) outage probability are derived for the case of perfect 
and imperfect CSI. These bounds are found to be reasonably tight and lead to bounds on the network 
transmission capacity. Third, the scaling laws for transmission capacity are derived for asymptotically 
small target outage probability. Specifically, with interference cancellation, the asymptotic transmission 
capacity grows according to a power law for both perfect and imperfect CSI. The base of the power law 
is the target outage probability, and the exponent is the inverse of the antenna-array size if it is smaller 
than the path-loss exponent. Otherwise, the exponent is bounded between the inverse of the antenna-array 
size and that of the path-loss exponent. Finally, the required lengths of CSI training sequences are derived 
for constraining the outage probability to be within a multiplicative factor of that with perfect CSI and 
for achieving a nearly optimal scaling law of the transmission capacity, respectively. These results are 
useful for evaluating the amount of CSI overhead needed to control the degree of network performance 
degradation due to CSI inaccuracy. 

Simulation results are also presented. As observed from these results, employment of a few (two to four) 
antennas per node is sufficient for harvesting most of the capacity gains promised by spatial interference 
cancellation. In particular, compared with the case of single antennas, a capacity gain of more than an 
order of magnitude can be achieved by using only one additional antenna at each node, even if CSI is 
imperfect. Moreover, a moderate length CSI training sequence is observed to be sufficient for keeping the 
loss of transmission capacity due to CSI inaccuracy small. These results demonstrate the effectiveness of 
spatial interference cancellation for practical applications. 

The remainder of this paper is organized as follows. Section |ll] describes the network and wireless 
channel models. Section |lll] introduces the algorithm for spatial interference cancellation, and presents the 
effective network and channel models based on the algorithm. The SIR outage probability and transmission 
capacity are analyzed for perfect and imperfect CSI in Sections ITVl and IVl respectively. Numerical results 
are presented in Section |VlJ followed by concluding remarks in Section I VII I 

II. Network and Channel Models 

A. Network Model 

In this paper, the locations of potential transmitting nodes in a MANET, including both active and 
inactive transmitters, are modeled as a Poisson point process following the common approach in the liter- 
ature [1], [2], [7], [36], [37]. Specifically, the positions of the potential transmitters form a homogeneous 
Poisson point process on a 2-dimensional plane with the density denoted by Aq. Time is slotted and 
in each time-slot potential transmitting nodes follow a simple ALOHA-like random access protocol, in 
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which they transmit independently with a fixed probability Pt. Let T„ denote the coordinate of the nth 
transmitting node. Given the random access protocol, the set ^ = {T„} is also a homogeneous Poisson 
point process but with the smaller density A = PjAq [38]. Each transmitting node is associated with a 
receiving node located at a fixed distance denoted as du 

Consider a typical receiving node located at the origin, denoted as Rq, and hence ITqI = d. This location 
constraint of Tq does not compromise the generality since the transmitting node process <l> is translation 
invariant. Furthermore, according to Slivnyak's theorem [39], the remaining transmitting nodes, namely 
$/{To}, remain as a homogeneous Poisson point process with the same node density A. 

The MANET is assumed to be interference limited and thus noise is neglected for simplicityo Con- 
sequently, the reliability of data packets received by the node Rq is determined by the SIR. Moreover, 
we assume that each data link in the network has a single stream, and communications between nodes 
are perfectly synchronized in symbols. All transmitting nodes are assumed to use uniform transmission 
power, Pf). Let S denote the random power factor for the link from Tq to Rq so that the received power 
at Rq is SPd- Likewise, let 5„ represent the interference power factor from transmitting node T„ to Rq, 
so that the received interference power at Rq due to T^ is In = SnPo- Thus, the SIR at Rq is given as 

SIR=^^ p. (1) 

^T„e<I'/{To} -''^ 
Since the SIR is independent of Pp, Pd = 1 is assumed for simplicity. The correct decoding of received 
data packets requires the SIR to exceed a threshold 6, which is identical for all receiving nodes. In other 
words, the rate of information sent from a transmitter to a receiver is no less than log2(l + 0) assuming 
Gaussian signaling. To support this information rate with high probability, the outage probability that 
SIR is below 9 must be no greater than a given threshold < e < 1, i.e. 

Pout(A)=Pr(SIR<^)<e (2) 

where Pout (A) denotes the SIR outage probability as a function of A. Given an outage constraint e, Pout 
determines the transmission capacity, which is defined as [2] 

C{e) = (1 - e)A, (3) 

where Pontile) = £■ Note that this equality maximizes transmission capacity under the outage constraint 
^ since Poutl-^e) increases monotonically with A^. 

^More generally, the distances between receivers and their corresponding transmitters could also be randomly distributed. 
However, as shown in [1], randomness in these distances has no significant effect on the analysis of transmission capacity and 
is thus omitted for simplicity. 

''Addressing the effect of noise requires straightforward but maybe tedious modifications of the analytical results in this work. 
In particular, accounting for noise changes only the linear factors in the asymptotic scaling laws for transmission capacity in 
Theorem |2] 
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B. Channel Model 

We adopt a narrow-band channel model with frequency-flat block fading, which is common in the 
literature (see, e.g., [2], [37], [40], [41]). During each transmission attempt the fading realization is 
assumed to stay fixed. Each node in the network is equipped with L antennas. Consequently, there 
exists a. Lx L multiple-input-multiple-output (MIMO) channel between every pair of nodes. Each MIMO 
channel consists of path-loss and spatially i.i.d. small fading components, corresponding to rich scattering. 
Specifically, the channel from a node r„ to the typical receiving node Rq is H„ = r„" G„. The factor 
r„" represents path-loss, where r„ = |r„| is the Euclidean distance and a > 2 is the path-loss 
exponent. The factor G„, is a L x L matrix of i.i.d. CM{0, 1) components, modeling spatially i.i.d. 
Rayleigh fading. Each transmitter is assumed to send a single spatial data stream so that stream control 
[42] is unnecessary and beamforming is applied at each transmitter and receiver. Let f„ and vq denote 
the beamforming vectors at T„ and Rq, respectively. Then the effective channel power for the data link 
from To to Rq is S = IvgHofoP, and that for the interference link from T„ to Rq is S'„ = IvgH^f^p. 

III. Spatial Interference Cancellation: Algorithm and Modeling 

The algorithm of zero-forcing beamforming for spatial interference cancellation is first described. The 
resultant effective network and channel models are discussed for perfect and imperfect CSI, respectively. 

A. Perfect CSI 

In this section, each receiver is assumed to have perfect CSI of the channel between each of its L 
strongest interferers and itself. 

1) Spatial Interference Cancellation and Opportunistic Transmission: The idea of spatial interference 
cancellation is to apply zero-forcing beamforming at Rq for canceling interference from strong interferers. 
The details of the interference cancellation algorithm are provided as follows. 

Let in and vq denote the transmit beamformer at T„ and the receive beamformer at Rq, respectively. 
From the perspective of Rq, the interference channel from T„ (n / 0) appears as an effective channel 
vector h„ = H„f„, where H„ denotes the actual MIMO channel. To facilitate our discussion, the indices 
of the transmitting nodes interfering with Rq are sorted according to their effective interference channel 
norms, namely ||hi|| > ||h2|| > • • • > ||h/,|| • • • . The crux of the interference cancellation algorithm is 
to constrain the beamforming vector vq of Rq to be in the null space of the matrix [hi, h2, • • • , h^^i]. 
Thereby, the interference from (L — 1) strongest interferers to Rq is nuUeqj: Ivghij = |vQh2|--- = 
|vQhi_i| = 0, where f represents the complex conjugate and transpose matrix operation. Note that 

^With probability one, the matrix [hi, h2, ■ ■ ■ , iiL-i] has full rank. Thus, with L antennas, _Ro can cancel at most (L — 1) 
interferers. 



July 11, 2008 



perfect CSI estimation of hi, h2, • • • , h^-i by Rq is required to completely cancel the interference from 
(L — 1) strongest interferers. CSI estimation at each receiver uses pilot symbols broadcast by transmitters. 
The issue of CSI inaccuracy is addressed in Section IIII-BI 

An arbitrary transmit beamformer is applied at Tq, represented by fp. Without performing interference 
pre-cancellation, Tq need not acquire CSI of the channels linking Tq and the interfered receivers, which 
is difficult. Note that an attempt to perform maximum ratio transmission [43] causes iterative updating 
of beamforming vectors at all nodes and potential network instability. By such beamforming, multiple 
transmit antennas contribute no diversity gain, but they are needed for interference cancellation when the 
transmitter becomes a receiver. 

To avoid deep fading due to the lack of diversity gain, opportunistic transmission is applied. Conse- 
quently, transmission at each transmitter is turned on only if the channel gain S = IvgHofoP is above 
a threshold denoted by /5o = f3d~°', where /? is an equivalent threshold appUed on the fading gain. It 
follows that the random-access probability for each potential transmitter is Pt = Fi{S > /3d~°) (cf. 
Section Hi- Al l. This algorithm is also used in [1] for single-antenna ad hoc networks and similar concepts 
exist in optimal power control for fading channels [44], [45]. The threshold /3 should be sufficiently small 
so as not to cause excessive delay before transmission]^ 

2) Effective Channel and Network Models: The following effective channel and network models result 
from the application of the interference cancellation algorithm in the preceding section. 

With perfect interference cancellation, Rq receives interference only from the node^ {T„ | n > L). 
Recall that r„ and /„ denote respectively the distance between T„ and the origin, and the interference 
power from T„ to Rq. Based on the channel model in Section Hl-Bl for n> L, the effective interference 
power In = PD|vjH„fnp = r^°|vjGn^p. Let /)„ = |vjG„f„p, so that for n > L, /„ = r^°p„. 
Because both f^ and vq are independent of G^ and G„ is an i.i.d. complex Gaussian matrix, the 
random variable /?„ will have an exponential distribution with unit mean. 

The effective power of the data link from Tq to Rq is given by S" = jVpHofoP = d^°|vgGofoP with 
S > (id"" due to opportunistic transmission. Recall that d is the distance between Tq and Rq, (3d"°' 
is the transmission threshold, and Go is the fading component of the MIMO channel Hq. Because the 
beamformers vq and fp are independent of Go as discussed in the preceding section, the random variable 

^Implementing this opportunistic scheme requires some form of hand-shaking between the transmitter and receiver and limited 
coordination among neighboring nodes. In particular, note that the transmission decision depends on the chosen value of vo, 
which in turn depends on the transmission decisions of the neighboring transmitters. We leave the development of such a protocol 
for future work. 

^Note that the indices of the nodes {Tn} are sorted according to the power of their interference to _Ro 
^Note that the interference canceling beamformer vo is a function of {G„ | 1 < n < L — 1} but independent of {G„ | n > L}. 
Moreover, the transmit beamformer f,i is random and hence independent of any interference channel. 
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Types of node 

W Transmitting node 
"5^ Receiving node 
l^ Primary interferer 
O Secondary interferer 
^Canceied interferer 




(a) Perfect CS I 



(b) Imperfect CSI 



Fig. 1. Effective channel and network models resulting from interference cancellation with (a) perfect CSI or (b) imperfect 
CSI for antenna arrays of three elements. The distance in the figures is proportional to the effective channel power. The data 
and interference links are plotted by using solid and dashed lines, respectively. 



W = IvJCofoP has the exponential distribution with the following probability density function 

fw{w) = exp{-w)/Pt, w> (3 



(4) 



where Pt = exp(— /3). 

Last, the effective channel and network models resulting from interference cancellation is illustrated 
in Fig. dja). Note that the density of transmitting nodes is A = PtXo- Based on the above models, the 
SIR at Rq is given as 



(Perfect CSI) SIR 



d-'^W 



Z^n=L ^" Pn 



(5) 



B. Imperfect CSI 

I) CSI Estimation: As discussed in Section JIIJ interference cancellation at Rq requires the estimation 
of CSI on the effective channels vectors corresponding to the (L — 1) strongest interferers of Rq. The 
CSI estimation is facilitated by the transmission of pilot symbols from these interferers. Let M denote 
the length of the pilot sequence sent by each interferer. The (L — 1) pilot sequences form the columns 
of a (L — 1) X M matrix, represented by Q, where M > {L — 1). Following [46], Q is designed 
as a unitary matrix. A protocol for CSI estimation may consist of two phases, where a receiver first 
identifies the (L — 1) strongest interferers by estimating their interference power and then estimates CSI 
for canceling these interferers. Given the unitary constraint for Q, this protocol should allow a transmitter 
to send different CSI training sequences to different receiversO The detailed design of the CSI estimation 



'Removing the unitary constraint for Q simplifies the CSI estimation protocol as a transmitter can broadcast the same training 
sequence to all receivers. Nevertheless, allowing the columns of Q to be non-orthogonal potentially introduces coupling between 
the estimated CSI for different interferers and complicates CSI estimation as well as analysis. 
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protocol will be addressed in future work. 

The CSI training signal received at Rq, denoted as y, is given as 

y^ = VM[hi,h2,--- ,hi-i]Q + 5^Kx^ (6) 

n=L 

where \/M ensures the transmission power of each node that sends pilot symbols is P^ = 1, and the 
M X 1 vector x„ contains CJ\f{0, 1) data symbols transmitted by the nth transmitter. The summation 
term in ^ represents interference to the CSI estimation at Rq. The CSI, denoted as hi,h2, • • • ,h.L-i, 
is estimated using the least-square method, thus 

r 1 1 -j^ oo 

hi,h2,--- ,hL_i =^=y^Qt = [hi,h2,--- ,hL_i] + ^= Vh„i^ (7) 

L J VM y^ n=L 

where x^ = x^Q^^ is a 1 x (L - 1) i.i.d. C7\A(0, 1) vector. Note that the alternative MMSE method 
requires estimating the covariance of the aggregate interference from the transmitters T^, T/^+i, • • • . Such 
estimation is potentially inaccurate due to the presence of the strong interferers Ti,T2,--- ,Tl-i. It 

follows from ([7]) that 

1 °° 
h/ = h^ -='^hnXn,i, I = l,2,--- ,L-l. (8) 



M . 

n=L 

The estimated CSI is applied for computing the beamforming vector vg used at Rq under the zero- 
forcing constraint vq -L [hi, h2, • • • , h^-i]. Thus, using ([U), the residual interference at Rq after beam- 
forming, denoted as Ir, can be written as 

L— 1 ^ L—1 oo 



t=l 1=1 n=L 



where a„ = r" vjhn is CAA(0, 1) since vq and h„ are independent. Therefore, I^ is CAA(0, o"^) with 
the variance a^ given as 



L-l 
.2 



-r. 



1=1 



/ J ^n O-nXrijl 



(9) 



n=L 

2) Effective Channel and Network Models: As illustrated in Fig.[TIb), the effective channel and network 
models for interference cancellation with imperfect CSI is identical to that for the case of perfect CSI 
except for the additional residual interference from the nodes {T„ | 1 < n < L — 1}. For the present 
case, the SIR in ^ is modified as 

(Imperfect CSI) SIR = tl^ -^^ ■ (10) 

'^i + l^n=Lrn pn 

It is worth mentioning that the present model of residual interference is more accurate than that in [7] 
for SIC with imperfect CSI. In [7], the residual interference power is modeled as z X^„=i r~°' where the 
parameter < z < 1 controls the degree of CSI accuracy. The overhead and algorithm for CSI estimation 
considered in the present model are not accounted for in [7]. 
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IV. Outage Probability and Transmission Capacity: Perfect CSI 

This section focuses on the analysis of the outage probability for a SIR constraint and the network 
transmission capacity for prefect CSI. In particular, the scaling law of transmission capacity is derived 
for the regime of small outage probability. 

A. Auxiliary Results 

To facilitate analysis, the interfering nodes of Rq after perfect interference cancellation, namely {T„ | 
n > L}, are separated into the strongest interferer T^ and the remaining interferers {T„ | n > L + 1}, 
referred to respectively as the primary and the secondary interferers. For convenience, denote the random 
interference power from Tl as G = II- There are two reasons for the above separation of the interferers. 
First, considering the strongest interferer Tl alone yields a lower bound for the outage probability to be 
derived in the next section. Second, the separation of interferers provides a useful result that conditioning 
on G, the secondary interferers {T„ | n > L + 1} form a Poisson point process as shown shortly. 

The result stated above is obtained by using the Marking Theorem [38]. To apply this theorem, a 
marked point process is defined for the secondary interferers, where the mark of the node T„ is the 
corresponding interference power /„. Specifically, conditioning on the interference power of the primary 
interferer G = g, the desired marked point process is 

U{g) = {{Tn,In} I Tn G <^/{To},0 <In<9] (H) 

where <l> is the homogeneous Poisson point process modeling all active transmitters (cf. Section III-AI ). 
Note that conditioning on G = g, the marks {/„} are independent. Given this condition, the result in the 
following lemma directly follows from the Marking Theorem [38]. 

Lemma 1: The process Il{g) is a homogeneous Poisson point process on M^ x M^ with the average 
number of nodes given by 

oo g 
^{g) =2ttX / rp{r,dl)dr (12) 


where p{r, •) is the distribution function of In conditioned on dn = r. 

This result is useful for analyzing the aggregate interference from the secondary interferers to Rq. 
Conditioned on G = g, this interference is written as 

Iu{g)= Yl ^"- (13) 

(T„,7„)Gn(g) 

The process luid) in (IT3] ) is known as a shot noise process [38]. The probability density function of 
luig) is difficult to derive and has no closed-form expression except for some simple cases [1], [37]. 
Nevertheless, using the Marking theorem, the first and second moments of this process are characterized 
in the following lemma. 
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Lemma 2: After perfect interference cancellation, the interference at Rq has the following properties. 

(a) The interference power of the primary interferer G has the following probability density function 

X L \L 

faig) = -Y^a-'''-' exp [-uXg-') (14) 

where 5 := - and u := TrT{l + 5). 

(b) Conditioned on G = g, the mean and variance of the aggregate interference power from the 
secondary interferers are given by 

niuig)] = ^/-' (15) 

Var(/n(5)) = ^,9'~' ■ (16) 

2 — 

Proof: See Appendix lAl D 

B. Bounds on Outage Probability 

From (Is) and the separation of interferers in the preceding section, the outage probability Pout can be 
written as 

Pout (A) = Pr(SIR <e) = E [Pr(/n(ff) >w%b-g\G, W)] (17) 

where tp := 9^^d~°', W is the fading component of the data link power for Rq (cf. Section Illl-A.2b . G 
and Inid) are respectively the interference power of the primary and secondary interferers. The direct 
analysis of the exact outage probability by using (|T7] ) is infeasible due to the difficulty in deriving the 
distribution function of Inid)- Instead, we bound Pout> following the approaches in [1], [2], [7]. 
The expression of the outage probability in (fTTll can be rewritten as 

PoutW = Pr(G > Wip) + E [Pr(In(5) > wip - g \ G < Wtp)] Pr(G < Wtp). (18) 

Thus, a lower bound of Pout is given as 

Pout{X)>'PT{G>WiP). (19) 

This lower bound considers only the primary interference, and hence is nearly tight if the primary interferer 
Tl is the dominant source of interference. Next, an upper bound on the outage probability can be derived 
by applying the following Chebyshev's inequality on dTSl ) 

Var(/n(g)) 
{a-E[In{g)]y 

Based on (fTSl ). (fT9l ) and (l20l ). bounds on the outage probability are derived as shown in the following 
theorem. 

Theorem 1: For perfect CSI, the bounds on the outage probability are given as follows. 



Pr(/n(5) > a) < min <j ^\;^^^\^ 1 |> , V a > E [Inig)] ■ (20) 
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1) The lower bound is 

plw = ^^^^ — - — TiLjpri ' ^"H-^^ '^ -V''- 

2) Define the sets Pi = {{w,g) \ < g < go,w > P} andV2 = {{w,g) \ g> go,w > g + E[Iu{g)]} ■ 
Moreover, let go denote a constant that satisfies the equation go + E[/n(5o)] = P^^ where E[/n((7) 
is given in Lemma |2] The upper bound on the outage probability is 

(w^-g-^^Xgl-S) 



P„V(A) = / / min ( , , '-' su , 1-3,2 ' 1 fw{w)fG{g)dwdg + P„(A) (21) 



where 



P„„) . l(^;:!0 _ 1^ r .--' exp f-,-V - A*-'A.- - .A.- ] .. 



r(L) PtT(L) ,,„ 

and 7(-, •) denotes the incomplete Gamma function. 
Proof: See Appendix iBl D 

The above bounds on Pout do not provide simple closed-form expressions in terms of the node density A. 
The difficulty in deriving such closed-form expressions is mainly due to the existence of multiple random 
variables, namely W, G and /n(G), which jointly determine the outage probability. The tightness of the 
above bounds on Pout is evaluated using simulation in Section IVll 

C. Asymptotic Transmission Capacity 

In this section, the scaling law for transmission capacity is derived for small target outage probability 
(e — > 0) and perfect CSI. This scaling law also accurately characterizes transmission capacity in the 
non-asymptotic outage regime (up to 0.1) as shown by simulations in Section IVll 

Small target outage probability results in a network of sparse transmitting nodes (i.e. A -^ 0). For such 
a sparse network, the useful relationship between the outage probability and node density is given in the 
following lemma. 

Lemma 3: For perfect CSI and A ^ 0, the outage probability scales with A as follows. 

1) For L < a, 

K, < lini ^2!^ < ^2 (22) 

A^O A^ 



where ki = -^ — p r(L+i) — ^"'^ '^2 = ki 



l + r^L,^ 



SLt 2 r,-S 



2) For L > a, 



Ki < lim — —f — , lim — ; < Ks (23) 

- A-^o A^ A^o A" ~ 



where k. - 4<5V'-^^^°r(-l /3)r(L-a+l) 
Where K3 - (2-S)T(L) • 

Proof: See Appendix O D 
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Note that the ratio ^ decreases as L reduces. This suggests that the asymptotic bounds are tighter for 
smaller values of L. 

Using Lemma [3] and the definition of transmission capacity in ^, the main result of this section is 
obtained and summarized in the following theorem. 

Theorem 2: For perfect CSI and small target outage probability e ^ 0, the transmission capacity 
scales as 

1) For L < a, 

lim-^>l, lim-^<l (24) 

where ki and K2 are specified in Lemma \3\ 

2) For L > a, 

lini-^>l, lim-^<l (25) 

where k^ is given in Lemma \3\ 

The above theorem shows that as the target outage probability decreases, transmission capacity grows 

following the power law ae* where a and t are constants. For L > a, only bounds on the exponent t 

are known. We conjecture that for L > a, t = 1/L. In other words, the exact scaling law is within a 

multiplicative factor of the asymptotic lower bound in the above theorem, which is confirmed by simulation 

results (cf. Fig. IHb)). The derivation of the exact scaling law for L > a requires a tighter upper bound 

on outage probability than that based on Chebyshev's inequality in (l20l ). This may require analyzing the 

distribution function of the secondary interference power (cf. Section IIII-A.2I) . which, however, has no 

closed-form expression for the present case [47]. 

For L < a, the exponent of the transmission capacity power law ae* is shown in Theorem |2] to be 

i_ i_ 

t = 1/L, and a is bounded as ^2 ^ < a < k^ ^ . This power law indicates that the size of antenna 

array L determines the sensitivity of transmission capacity to the change on the outage constraint. To 

facilitate our discussion, rewrite the scaling law in Theorem [2] as C(e) = ae~ where "=" represents 

asymptotic equivalence for e ^ 0. Moreover, consider two sets of values (Ci,ei) and (C2,e2), and 

define the logarithmic ratios AC = log ^ and Ae = log |^. Using this notation, the above scaling law 

can be written as 

The above quantity ^ represents the sensitivity of transmission capacity towards changes in the outage 
constraint. Its value decreases inversely with the size of antenna array. Specifically, computed using 
(l26l ). a hundred-time decrease on e reduces network transmission capacity by {10,3.2,1.8} times for 
L = {2, 4, 8}, respectively. For the extreme case of L = oo, transmission capacity is independent of the 
outage constraint since -^ = 0. Last, from simulation results in Section |VlJ the capacity scaling law in 
Theorem |2] is observed to also hold in the outage regime of practical interest (e < 0.1). 



^out (A) = Pr ( , . _^_i _ < e ) (27) 
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V. Outage Probability and Transmission Capacity: Imperfect CSI 

In this section, the SIR outage probability and transmission capacity are analyzed for the case of 
imperfect CSI. 

From ( ITOl ), the SIR outage probability for the case of imperfect CSI is written as 

where a"^ is the power of the residual interference given in ^. This probabiUty is related to that for the 
case of perfect CSI in the following theorem. 

Theorem 3: For the case of imperfect CSI, the SIR outage probability is bounded as 

Pout(A)<Pout(M,A)<Pout(A)(l + f)+2^"ie-'^^, M>L-1 (28) 

1 
where Pout is the outage probability given perfect CSI, Z > is arbitrary and to = [T{L)] ^-i. 

Proof: See Appendix iDl D 

The result in Theorem [3] can be interpreted using two corollaries. The following corollary is obtained by 
setting Z = M« in (|28l) with q £ (0,1) and M ^ cx), and then let q ^ 0. 

Corollary 1: For M — > oo and fixed A, the outage probability given imperfect CSI converges to that 
for perfect CSI as follows 



Po„t(A) \M, 

This corollary agrees with the intuition that increasing the amount of CSI training overhead reduces the 



effect of CSI inaccuracy. Moreover, (I29b shows that the effect of imperfect CSI disappears at a rate 
proportional to the inverse of the pilot sequence length. The next corollary of Theorem |3] characterizes 
how much CSI training overhead is needed for containing the effect of CSI inaccuracy. 

Corollary 2: To ensure Pout{M, A) < (1 + ^)Pout(A) with ^ > 0, it is sufficient to choose the length 
of the CSI training sequence as 



M = max 



^(Llog2-logPo,t(A)-loge) 



L-1] . (30) 



where \a\ gives the smallest integer larger than a. 

This corollary is derived by substituting Z = —^log (2^^PoutS,) and M = 2Z/^ into (1281) . and taking 
into account that M > L — 1 and M is an integer. The result in Corollary [2] shows that for given ^, M 
increases linearly with L and logarithmically with the inverse of Pout- 

Finally, the required length of pilot sequence is derived for achieving the scaling laws of transmission 
capacity close to those for the case of perfect CSI. Define the transmission capacity for imperfect CSI 
as C{M, e) = (1 — e)P^J^{M, e). The scaling laws of C{M, e) are given in the following theorem. 
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Theorem 4: For < (/? < 1, choosing the length of pilot sequence for CSI estimation as M^{lp) 



max 



M{^) 



, L - 1 1 with 



M{ip) 



"^ [(L-l)log2-21oge], L<a 



-— ^^— [(a-l)log2-21oge], L>a 

achieves the scaling laws of transmission capacity for a small target outage probability (e -^ 0) given as 
1) For L<a: 



2) For L > a: 






e— -^ J- ~ e^O -r J- ~ 



where ki, K2> and k^ are given in Lemma [3l 

Proof: See Appendix El D 

Theorem ^ shows that CSI inaccuracy introduces only an additional linear throughput scaling factor, 
namely < 93 < 1, without affecting the exponents in the throughput scaling laws. As shown by Fig. [2l 
the required length of pilot sequence M^ increases relatively gradually with growing ip in the regime of 
< (/? < 0.6 but very rapidly in the regime 0.6 < (/? < 1. In the regime of small ip, M^ grows as L 
increases due to the constraint that M > L—1. Nevertheless, the reverse holds in the large tp regime. The 
reason is that increasing L reduces the number of strong interferers for CSI estimation at each receiver, 
and hence requires less CSI overhead. 

VI. Simulation and Discussion 

In this section, the bounds on outage probability and the network transmission capacity are evaluated 
using Monte Carlo simulation. The procedure for simulating a MANET follows that in [48]. The simulated 
ad hoc network lies on a two-dimensional disk and contains a number of transmitter-receiver pairs, which 
follows the Poisson distribution with the mean equal to 200. The locations of the nodes are uniformly 
distributed on the disk. The disk area is adjusted according to the node density. The typical receiver is 
placed at the center of the disk. We set the distance between the typical transmitter and receiver as d = 5 
m, the required SIR as = 3 or 4.8 dB, the transmission threshold /3 = 1, and the path-loss exponent as 
a = 4. 

A. Bounds on Outage Probability 

For perfect CSI, the bounds on outage probability from Theorem [J and simulated values are compared 
in Fig. [3] The number of antenna per node is L = {2, 4}. As observed from Fig. [3l the bounds for L = 2 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Linear throughput scaling factor, (t> 



Fig. 2. For imperfect CSI, the required length of pilot sequence Afj vs. the increasing throughput scaling factor if for 
L = {2,3,4} and e = 0.01. 

are tighter than those for L = 4. Moreover, the bounds and the simulated values of the outage probability 
converge as the transmitting node density A decreases. These two observations can be explained by the 
dominance of the primary interference over the secondary one as L increases or A decreases, where 
the secondary interference causes the looseness of the bounds on outage probability. Finally, the outage 
probability is approximately proportional to A^. 

In Fig. m the SIR outage probability for imperfect CSI is evaluated against that for perfect CSI, where 
the transmitting node density is varied and different lengths of pilot sequence are considered. The number 
of antennas per node is L = 4. The length of pilot sequence is fixed at M = {3, 5, 11} in Fig. IHa) but 
varied with the node density in Fig. |4jb) based on Corollary |2] As observed from Fig. Ufa), increasing M 
rapidly converges the outage probability for imperfect CSI to its lower bound corresponding to perfect 
CSI. In particular, for M = 11, CSI inaccuracy increases outage probability by less than a factor of two. 
In Fig. HJa), M varies with the node density according to (|30l ) with ^ = {1,5}. Regardless of the node 
density, the outage probability for imperfect CSI is bounded within (1 + C) times its lower bound, which 
validates Corollary |2] For instance, the ratio between the outage probabilities for imperfect and perfect 
CSI is about 3 with is smaller then the upper-bound (1 + ^) = 5 given in Corollary |2] 



B. Scaling Laws of Transmission Capacity 

In Fig. [51 asymptotic bounds on transmission capacity in Theorem |2] are compared with the exact 
values obtained by simulation for perfect CSI and the range of target outage probability e G [10^^,0.1]. 
The corresponding curves are identified using the legends "asymptotic upper bound", "asymptotic lower 
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10 10 

Transmitting node density 



(a) Perfect CSI, i = 2 



Transmitting Node Density 
(b) Perfect CSI, i = 4 



Fig. 3. Outage probability for different transmitting node densities and perfect CSI. The size of the antenna array is (a) L = 2 
and (b) L = 4. 

bound", and "simulation". Different combinations of L = {2, 4, 5} and a = {3, 4} are separated according 
the cases of L < a and L > a, corresponding to Fig. [5la) and Fig. Ob), respectively. As observed from 
Fig. Oa), for L < a, the asymptotic upper bound on transmission capacity is tight even in the non- 
asymptotic range e.g. e G [0.01,0.1]. The tightness of this bound is due to the dominance of primary 
interference for interference cancellation with small sizes of antenna array. Moreover, Fig. |5lb) shows 
that for L > a the slopes of the "simulation" curves converge to those of the corresponding "asymptotic 
upper bound" curves as the target outage probability decreases. The above observations suggest that for 
both L < a and L > a, the scaling laws of transmission capacity for small target outage probabilities 
follow the power laws with the same exponent equal to l/L. 



C. Transmission Capacity vs. Size of Antenna Array 

In Fig. |6j the transmission capacity is plotted for an increasing number of antennas per node assuming 
perfect CSI. Furthermore, different outage constraints, namely e = {10~^, 10~^, 10~^}, are considered. 
From Fig. |6l the following observations are made. First, the use of multiple antennas for interference 
cancellation leads to the increase in transmission capacity by an order of magnitude or more with respect 
to the case of single-antenna per node. This capacity gain is especially large for a small number of 
antennas and small target outage probability. For example, for e = 10~^, the use of three antennas per 
node provides transmission capacity seven times of that for the single-antenna case. The capacity gain 
by using additional antennas diminishes rapidly as the number of antennas per node increases. Second, 
the outage constraint affects transmission capacity significantly for a small number of antennas per node. 
Nevertheless, transmission capacity becomes insensitive to changes in the outage constraint as the number 
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Transmitting node density 



Transmitting node density 



(a) Fixed M (b) Adaptive M 

Fig. 4. Compare outage probability for perfect and imperfect CSI given different transmitting node densities. The length of 
pilot sequence is (a) fixed at M = {3, 5, 11} or (b) adapted to the node density according to J30b for ^ = {1, 5}. The size of 
the antenna array is L = 4. 

of antennas increases. 

The effect of imperfect CSI on transmission capacity is shown in Fig. |71 where transmission capacity 
is plotted for different numbers of antennas per node L. The size of the antenna array is L = 4 and the 
outage constraint is e = 10~^. The length of pilot sequence is M = {3, 5, 11}. Several observations can 
be made from Fig.|7] First, the loss on transmission capacity due to CSI reduces as M increases. Such loss 
is relatively small even for a moderate value of M. For instance, the reduction on transmission capacity is 
25% for M = 11 and L = 8. Second, even for minimum CSI estimation overhead (i.e. M = 3), a capacity 
gain of more than an order of magnitude can be achieved using interference cancellation. This supports 
the practical applications of interference cancellation. Finally, most capacity gains are contributed by the 
cancellation of the strongest interferer to each receiving node. The cancellation of more interferers has a 
much less significant effect on the network capacity since it becomes limited by residual interference. 



VII. Conclusion 

In this paper, a spatial interference cancellation algorithm is applied in a MANET setting and the 
resulting gains in the network's transmission capacity are characterized under given constraints on the 
SIR outage probability. Bounds on this outage probability are given for a network with Poisson distributed 
transmitters. For asymptotically small outage probability, the scaling laws of transmission capacity are 
derived, which follow the power law for both perfect and imperfect CSI. These scaling laws also 
accurately predict transmission capacity for the non-asymptotic outage regime. The required lengths of 
CSI training sequence are derived for constraining the increase of outage probability due to CSI accuracy, 
and for achieving close-to-optimal capacity scaling laws, respectively. Through simulation, interference 
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Fig. 5. Comparison between asymptotic bounds on transmission capacity and the exact values obtained by simulation for 
perfect CSI and the cases of (a) L < a and (b) i > a. 

cancellation are observed to provide significant network capacity gains even for just a few (2 — 3) antennas 
per node. As also observed, a moderate amount of CSI training overhead can effectively alleviate the 
impact of CSI inaccuracy on network transmission capacity. Thus, spatial interference cancellation may 
be suitable for appUcations in practical network. 

This work opens several issues for future work. First, the bounds on the SIR outage probability derived 
in this paper provide no simple closed-form expressions. Furthermore, characterizing the asymptotic 
behavior of the outage probability also seems difficult. Thus, new and more powerful analytical tools 
should be found if possible. Second, integration of interference cancellation and stream control [42] 
provides additional gains on network transmission capacity. Last, it is useful to investigate the performance 
of interference cancellation for ad hoc networks in more realistic settings such as heterogeneous traffic 
patterns and high mobility. 
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Fig. 6. Transmission capacity by simulation for different node densities and perfect CSI. The size of the antenna array is L = 4 
and the outage constraint is e = {10^^, 10^^, 10~'^}. 
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Fig. 7. Transmission capacity by simulation for different node densities and imperfect CSI. The size of the antenna array is 
L = 4 and the outage constraint is e = 10~^. 
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Appendix 

A. Proof of Lemma |2] 

Define a marked Poisson point process as M-{g) = {(r„,I„} | r„ G <i'/{To},I„ > g}. Given that 
In = r~'^pn, the node density of M.{g) follows from the Marking Theorem as 

lx{M{g)) = X / 27Trfp{t)drdt 

Jo Jo 

/•oo 

= nXg-^ / t^e-^dt 



JO 

= yXg-^. (32) 

where 5 and v are defined in the statement of the lemma. The cumulative density function of G is the 
probability that the number of nodes in the subset M^{g) is no more than (L — 1). Thus, 

L-l / ■, -5\fc 

Pr(G <g) = Y. ^! ^""P (-^^5"') • (33) 

fc=0 

The probability density function of G is obtained by differentiating the above function 



^ \v\f+\_^,^^-^.-^^_^y^„-s ^"^-A {vxf 



foia) = SY^'-^^g-'^'^'^-'e-^'^-'-SY^ 



Sk-l -u\g 



k=0 k=l ^ ' 

lfc=o ■ fc=o ■ J 

The desired result follows from the last equation. 

Using Campbell's theorem [39], the expressions for Var(/n((7)) and E[In((7)] are obtained as 



E[InU] = 2^A / \ y~^pe-Pdrdp 
Jo Jin 



(f)' 



^9-'. (34) 

a — 2 



Var(/nU) = 2^A / l^.r{r"'p)\'Pdpdr 
Jo J{^) 

f'OO f'OO 



^g'-' r p'e-^dp 
a- I Jo 



"' .-^ (35) 



a — 1 
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B. Proof of Theorem |7] 

Define the product space of {W, G) as V = [/?, oo) x [0, oo), and ip = 9~^d~°. The outage probabiUty 
can be written as 

^out = J J Pr(In(5) +5 > w)fwHfGig)dwdg. (36) 

As illustrated in Fig. [8j "D is partitioned as D = U^^qD„ where Di and V2 are defined in the statement 
of the theorem and 

2?o = {{w,g)\0<g>wi>,w>P} (37) 

^3 = {{w,g)\go<g<Pil^,p<w<g + B[In{g)]} (38) 

^4 = {{w,g)\g>P^,g<w<g + -E[In{g)]}. (39) 

In Vq, the primary interference G is sufficiently large for causing an outage. Considering Vq alone gives 
an lower bound on the outage probability denoted as P^^^^ and given as 

^out= // PT{In{g)+g>w)fw{w)fG{9)dwdg= jj fw{w)fG{g)dwdg. (40) 

(w,g)&Vo {■w,g)e'Do 

The subset Vi U P2 satisfies the criterion for applying Chebyshev's inequality in (|20l ). In D3 U V4 where 
Chebyshev's inequality does not hold, the outage probability is upper bounded as Pr(/n(5') +g > w) < I. 
Based on the above discussion and (l20l ). an upper bound on the outage probability, denoted as P^ut^ i^ 
written as 



{■w,g)eViUV2 



(41) 



n 
fwiw)fG{g)dwdg. 

{w,g)£V3UV4 

The bounds in (l40l ) and (|4T] ) are further developed in the following sub-sections. 
1) Lower Bound on Outage Probability: By substituting ^ and (IT4l ) into (l40l ) 

Si,L\L coo i-glii) 
PL = ^ / e-g-'^-'e~^>^^-'dwdg (42) 

s,,,L \L r /■CO 1 />oo 

= ^^ / g-SL-l^-.Xg-^^g_^ g-'^-'e-'^^^a-^-^-'^dg 

r(L) r(L)p, y, ' ""'■ 

The last equation gives the desired lower bound on the outage probability. 
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Fig. 8. Partitions in the product space V (shaded) of (W, G), which are used to obtain bounds on the outage probability. 



2) Upper Bound on Outage Probability: From (|4T 



^o"ut = PLt+^+ I fwHfG{9)dwdg 

n+ I fw{w)fG{g)dwdg 



VL+ I I fw{w)fG{g)dwdg 

Jgo Jf3 



Q + 



du^X^ rg-SL-l^-uX,-^^g_^J^ rg-SL-l^^i'-^9-^^^-^X9^-^-^X9-^dg_(^^^^ 



r(L) 



9o 



^{L)Pt 



ga 



The desired upper bound follows from the last equation. 



C. Proof of Lemma\3\ 

1) Asymptotic Lower Bound on Outage Probability: For A — > 0, with cq defined following (IT4l ). the 
lower bound on the outage probability in (l42l ) can be rewritten as 

PL. = "„'1V:. / / g'^^'^ exp (-i^Xg-^) e-'"dgdw 



Pt^{L) J p Jwip ^ ' 

^-^r r9-''-'[l + 0iX)]e-dgdw 

PtT{L + I) jp 

KiA^ + 0(A^+^), A^O 



w-^^e-'^dwX^ + 0{X^+^), A^O 



(44) 



where ki is given in the statement of the lemma. 
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2) Asymptotic Upper Bound on Outage Probability: Define a set and its complementary set, represented 
by 2?5 and Dg, as follows 

P5 = JKS) I W > fj^g + -^—^g^-^ < -^\ , VI = Uw,g) I W > (i,g+-^-—^g^-^ > -^ 

Moreover, let g{w) and g{w) denote the functions that satisfy g + j^A^^^*^ = ^ and g + j^A^^"'' = 
wil), respectively. Note that "Ds G Pi U I?2> where "Di and 1)2 are defined in Theorem [1] By using the 
above definitions, the term Q. defined in (|4T| | can be upper bounded as 

Var(/^(g)) 

,(u;V-<7-E[/.(<7)])^ 



^ < // ^^^[ ....,,. ' "^r. .,x.,2 >l)/H-M/G(g)^^^g + 



{»i;,3)e{X'iUX'2)nD5 

fw{w)fG{g)dwdg 






/3 JO 
oo /•g(w)) 

/ fw{w)fG{g)dwdg 

(a) /-oo |.3(«,) 4Jj^;^^2-5 roc .g(w) 

< / / 77, — K7 — 7^fG{9)fw{w)dgdw+ / fw{w)fG{g)dwdg, A ^ 0(45) 



■• ^_ 



A2(A) Ai(A) 

Obtaining the inequality (a) uses Lemma |2] and the definition of g. By observing that g < wip/2 from 
its definition and substituting (fT4l) . A2(A) defined in (1451 ) is upper bounded as 

A2(A) < ^^^/°°u;-^^'^^<7-^(^+^)+^exp(-.A,-^)d,/H.Hdu; 

= 7^ — JyF77\ / ^"^ / g^~" eyiv{-g)dg fw{w)dw. (46) 

(2 - (^)r(L) 7^ Ju\{^)-^ 

'■ • V ' 

A3 (A) 

Note that for A ^ 0, A3 (A) in (l46l) is bounded for L > q and unbounded for L < a. It follows from 

(l46l) that 

r(L-Q + l), L> a 

-5L+2 (47) 



xL-Q / w^p 



A3(A) < <( 

By substituting (07]) into (l46]l 

(2 - 5)r(L) 



(^A)"-(_l) [1 + 0(A)], L<a. 



' 45^-2i/°E [W-^] r (L - a + 1) ^ 

A , L/ ^ OL 



A2(A) < <^ 
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Next, an asymptotic upper bound on Ai(A) defined in ( |45] ) for A — > is obtained as follows. Since 

v-ip 

T 



g{w) < wtl) and g{w) = ^ + 0(A) from their definitions, by substituting ((14)) . Ai(A) is bounded as 



Ai(A) < 



6u^X^ roo rw4^ 



9 



'SL-1 



< 



2 



exp{—i'Xg )dgfw{w)dw 



^L-l 



g e^p{-g)dgfw{w)dw 



j,A(^)-5 + 0(A2)](i-i) 



/3 



r(L) 



exp(— i/A(tt;V')' 



exp(-z.A(^)-^ + 0(A^)) 



fw{w)dw 



i,i-iA^-i(!^)-^(^-i)+0(A^) 



-S/nS 



uX{w^y'{2 



1)+ 



r(L) 



+ 0(A 



L + lN 



By combining 



48] ). and ( |49l ), an asymptotic upper bound for Q, in (|4TI ) is obtained as 

' 45V'"^z^"E [Ty-2] r (L - a + 1) 



J]< < 



(2 - 5)r(L) 



-A", 



L > a 



^5L 



^V'^^E [T^- 



5L1 



-A^ + 0(A^+i), L<a 



(49) 



(50) 



,2 -(5 J r(L) 

where E[W^-'5^] = r(l - 5L,fi)/Pt and E[VF-2] ^ r(-l,/?)/Pj. 

Finally, using (70 = /'V' + 0(A) and following the similar procedure as in Section ICTl it can be shown 
that il(A) = KiA^ + 0(A^+^), where Vt is defined in (|43]). Combining this result, (|45]) and dSO]! leads to 
the desired asymptotic upper bound on the outage probability. 



D. Proof of Theorem \3\ 

Using (ITOl ) and by definition, the conditional SIR outage probability is given as 



Pr(SIR < 6) 



E 
E 



Pr(SIR<0|$,{/5„},W^) 



Pr ( Wd-'^O'^ <(tI + Y^ 



^n Pr 



n=L 

00 



^,{Pn},W 



L-l 



= E PrI Tyd-"0-i<l:^r-V„^z, + ^r->„ cI,,{p4,vfJ (51) 

\ n=L 1=1 n=L / . 

where the i.i.d. random variables {zi} follow the exponential distribution with unit mean. Let (^ denote a 

chi-squared random variable having (L — 1) complex degrees of freedom. It follows from (ISTI ) that 

00 



Pr(SIR < 61) = E 



Pr(wd-^9-'<f2r-yJl + j^ 

\ n=L ^ ^ 



'^,{Pn},W 



(52) 
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The above equation can be expanded as 
Pr(SIR <e)= B 






'^,{pn},WX>Z\ 

<^,{Pn},W,C<Z 



Pr(C > Z)+ 



Pr(C < Z) 

(53) 

Using the key result in [49], the incomplete Gamma function with (L — 1) complex degrees of freedom 
can be upper bounded as 



r(L - 1) Jz 

where u; = [T{L)]'^^^. By combining ^3^ and (|54] 
Pr(SIR<0) < E 



T^-^e-^dT < 1 - (1 - e""^^) 



,L-1 



E 






<l>,{Pn},W,C<Z 



eC^^'h-^y- 



n=L 

l-LUlZ 



<^,{pn},W,C>Z 



< E 



Pr Wd-'^e-' < Y, r-^n 1 + V7 



n=L 



M J 



^APn},W 



+ 



^Y^T^U-lVe--'^ 



1=1 



< E 



FT[Wd--9-'<Y^r-yJl + jj) 



n=L 



z\ 



M 



^APn}.W 



+ 



EP/V-iye--'^ 



[=\ 



I 



Since W follows the exponential distribution. 



Pr(T^ < a(l + h))= / e'^dr < (1 + 6) / e-^dr. 

Jo Jo 



From (|55]l and (|56]l 

Pr(SIR < 61) < E 



Pr Wd-'^e-^ < f2 ^n Vn 



n=L 



'^,{Pn},W 



(54) 



+ 



(55) 



(56) 






E f^ 7 'V-D'^-""^ 



«=1 
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L-1 



H-ujIZ 



pr(siR <«) (i + 1) +.-^X;(^ ;')(-!)' 

Pr(SIR <e)(l + — \+ 2^~^e''^^. (57) 



The desired result follows from dSTl) . 

£. Proof of Theorem E] 

By substituting Z = -^^ log (21"-^ e^) and M, = 4 into (EH) 

Pout(M„A)<Pout(A)(l + (^) + e^ (58) 

To simplify notation, denote Pout {Me, A) as Pout (A). Define A^, A^ and A^ using the equations Pout (K) = 
e> ^out(A") = e and Pout(Ae) (1 + <^) = e — e^> respectively. These quantities satisfy the relationship 
a' < A < A". The first inequality is due to that the functions Pout (A) and Pout (A) monotonically decrease 
with increasing A; the second holds since CSI inaccuracy reduces the maximum transmitter density under 
the outage constraint. For L < a, using the above inequalities and Theorem |2] 

Al 



1 _£(lz£) 



> (1 + (^)-i/^ (59) 



lira ^ — > lira 

,i„.£«ll) < ,i„(l^<l, (60) 



By substituting y? = (1 + ipyT into ([591), ([601) and M^ = - J^ log (2^~^e^) and taking into account that 
Me > L — 1 is an int 
the same procedure. 



ipu> 

Me > L — 1 is an integer, the desired result for L > a follows. The result for L > a is derived following 
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